Application of indirect Hamiltonian tomography to systems with short coherent times 
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The identification of parameters in tiie Hamiltonian that describes complex many-body quantum systems 
is generally a very hard task. Recent attention has focused on such problems of Hamiltonian tomography 
for networks constructed with two-level systems. For open quantum systems, the fact that injected signals 
are likely to decay before they accumulate sufficient information for parameter estimation poses additional 
challenges. In this paper, we consider use of the gateway approach to Hamiltonian tomography |T|,[2|] to complex 
quantum systems with a limited set of state preparation and measurement probes. We classify graph properties 
of networks for which the Hamiltonian may be estimated under equivalent conditions on state preparation and 
measurement. We then examine the extent to which the gateway approach may be applied to estimation of 
Hamiltonian parameters for network graphs with non-trivial topologies mimicking biomolecular systems. 



I. INTRODUCTION 

Precise information about the Hamiltonian of many-body 
quantum systems is crucially important for analysis and pre- 
diction of their dynamics, especially to understand the extent 
to which a given subsystem behaves quantum mechanically. 
If the subsystem of interest is well isolated from the remain- 
der, i.e., from its environment, in the sense that its dynamics is 
immune to the effect of noise, then the time evolution is uni- 
tary and observable data may be considered 'clean' enough 
to extract good information on the subsystem Hamiltonian. 
However this is clearly an idealized situation that is rarely en- 
countered. Furthermore, the general procedures to estimate 
even a small part of the Hamiltonian that acts on a many-body 
system are, in general, very complex and require a large num- 
ber of measurements of different observables. This leads to a 
further challenge which is that as the data acquisition process 
becomes more elaborate and accesses more of the system, it 
is likely to introduce an increasing amount of measurement 
noise. 

Several approaches for Hamiltonian identification, or 
Hamiltonian tomography, have been recently proposed that 
seek to reduce the complexity of the procedure by making 
use of some a priori knowledge about the physical system yl- 
[Tlll . One approach is to map the many-body system onto a 
quantum network and to make use of knowledge about the 
topology of this network to devise protocols that extract de- 
sired Hamiltonian parameters from measurements on a re- 
stricted portion of the network. Following the demonstration 
that the Hamiltonian parameters of one-dimensional chain of 
spin- 1/2 particles may be determined by measurements on 
a single spin iH this approach has been generalized to 
more general spin networks with restricted measurement ac- 
cess on a small gateway region fi\ as well as to more general 
Hamiltonians lfl3ll . Reference 112J also showed that such an 
estimation scheme may be robust against noise under weak- 
coupling conditions. For sparse Hamiltonians, a different ap- 
proach has recently been developed using the method of com- 
pressed sensing llT4[[l5[ ]. which has also been applied to quan- 
tum state tomographv J l6ll . The compressed sensing approach 



allows determination of both higher order Hamiltonians and 
system-bath interactions, but is limited to sparse Hamiltoni- 
ans. Other approaches have been developed based on convex 
optimization iITtIi and Bayesian estimation i^. While these 
Hamiltonian tomography approaches are related to the bet- 
ter known quantum process tomography (QPT) lfl"8l - l22ll that 
(together with quantum state and quantum measurement to- 
mographies) provides a complete characterization of quantum 
dynamics, they differ from QPT in seeking to reconstruct the 
desired parameters with a minimal amount of resources. 

In this paper we explore the use of the gateway scheme 
outlined in 121] for determination of Hamiltonian parameters 
for an open quantum system under conditions of restricted ac- 
cess. The approach of Ref. |2|] was based on the assumption 
of long coherence times, which allowed the injected signal 
(spin wave) to go back and forth in the network many times so 
that the information about spin interactions may be encoded in 
the signal. For dissipative systems, we cannot in general ex- 
pect such a long lifetime of the signal and it will generally be 
susceptible to decay before coming back to the injection site, 
even though the initial time evolution for a short time may be 
seen to be coherent. We therefore limit our attention here to 
complex systems in which a subsystem does show such co- 
herent short time evolution. 

One prototype of this latter situation that is of consider- 
able cuiTent interest is the subsystem of pigments in photo- 
synthetic light harvesting systems. Recent spectroscopic ex- 
periments have shown that electronic energy transfer dynam- 
ics in such systems displays coherence for several hundreds of 
femtoseconds ll23l426ll . Although effects of vibrational contri- 
butions are not entirely clear, these coherences are generally 
accepted to reflect quantum coherences between different ex- 
citonic states that may be described by superpositions of sin- 
gle molecule electronic excitations, and are thus amenable to 
a two-level pseudo-spin representation. In this work we shall 
consider a network of pseudo-spins that mimics pigments in a 
light harvesting protein. 

We first review the gateway scheme of Refs. [Ql, 01, in- 
troducing the graph theoretic description and notion of infec- 
tion between different regions of the network (graph) (Sec.HIb. 
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We then summarize the minimal restrictions on measurement 
access via spectroscopic measurements in a pigment-protein 
complex (Sec. HiH i. These differ from the measurement re- 
quirements for spin networks and thus necessitate an ex- 
tension of that approach. We present a classification of net- 
work topologies that are accessible to Hamiltonian tomogra- 
phy under the current scheme. In Sec.|IV]we then investigate 
the extent to which the scheme may be applied to a network 
graph mimicking pigments embedded in the Fenna-Matthews- 
Olson protein of photosynthetic green bacteria iIztIi . A discus- 
sion and analysis of the limitation posed by restriction to short 
time scales of coherent evolution, together with indications for 
extensions to remedy this, follows in Sec. [V] 



II. GATEWAY SCHEME OF HAMILTONIAN 
TOMOGRAPHY 

In the 'gateway scheme' of Hamiltonian identification 
JUS], we consider a network of spin 1/2 pseudo-spins sub- 
ject to a unitary dynamics generated by a Hamiltonian con- 
taining pairwise interaction terms and Zeeman terms. For 
clarity and conciseness we describe here only the case of 
excitation-conserving Hamiltonians, namely those satisfying 
[H, Zn] = 0, i.e., conserving the total magnetization of 
the pseudo-spin network. We further assume that all coupling 
strengths c„ between spins are real and the (relative) signs of 
these, but not the magnitudes, are known. We shall consider 
the determination of the Hamiltonian parameters in the first 
excitation subspace, i.e., the subspace in which there is only 
one 'up' pseudo-spin and all other pseudo-spins are 'down'. 
(Note that this places a restriction on the interaction between 
pseudo-spins.) We illustrate the scheme here for a ID spin 
chain with nearest-neighbor interactions, for which the Hamil- 
tonian in this subspace is given by 
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For a ID chain, we start our procedure by measuring 
all eigenenergies {Ej} and coefficients {(£'j|l)} for j = 
1, 2, . . . , A^. Here \Ej) are the energy eigenstates, |1) denotes 
Itl ... I) = |10 . . . 0), and N is the number of pseudo- 
spins. Other states with a single up-spin will be denoted 
similarly hereafter, i.e., |n) = || ... ^ • ■ • i) contains 
only a single spin up at the n-th site. We define Pn{t) = 
tr^„ [[/(<) I V'o) {'>Po\U{ty] to be the reduced density matrix on 
site 71, where U (t) = exp{—iHt) is a time evolution operator 
for the chain. 

First we initialize the state of the chain to be |-0o) ■= 
l/v^(|0) + |1)), namely, l/\/2( 1 0) + |0 . . . 0). Such an 
initialization is possible by accessing only the first spin ll28ll . 
We then perform state tomography on the first spin after a 
time lapse t to extract the reduced density matrix pi (t) and 
repeat this at various time delays to obtain pi (t) as a function 
of time. Up to an irrelevant phase factor, the diagonal element 



of pi (t) can be written as 
fnit) := (l|expHiJt)|l) = ;^cxp(-zS,t)|(^j|l) 



The eigenenergies {Ej} and the coefficients can 
then be obtained by performing the time Fourier transform 
of fn{t). Due the arbitrariness of the global phase, we can 
choose all {Ej\l) to be real and positive. Detailed discussion 
of the factors determining the efficiency of the Fourier trans- 
form are discussed in Ref. UJl. As noted there, it is neces- 
sary to observe the repeated reflections of the signal (at least 
N times) in order to obtain an adequate signal to noise ra- 
tio. Consequently a long time coherence is necessary for im- 
plementation of Hamiltonian tomography with measurements 
only on a single spin. 

With the information about {Ej} and {{Ej\l)} obtained 
from these single spin measurements, we can proceed the pa- 
rameter estimation by constructing a set of N'^ equations rep- 
resenting {Ej\Hiu\n) for 1 < j,n < N. These equations are 
given by 



E,{E,\l)=h{E,\l)+c,{E,\2), 
Ej{Ej\n) = Cn-i{Ej\n - 1) + 5„(i?,|n) 

+ Cn{E,\n+l) (KjkTV), 
E,{E,\N) = CN-i{E,\N - l)+bN{E,\N). 



(2) 

(3) 
(4) 



Noting that bi ^ = each factor of 

which has been determined by state tomography on spin 1 , the 
expansion of |2) in the basis \Ej) can be obtained up to the 
constant c\. The value of c\ can then be obtained within a sign 
factor by requiring that |2) be normalized to unity. All other 
parameters 6;, Cm and {Ej\ri) are then subsequently obtained 
in the same manner from Eq. (O (for I = N —l,m = N) and 
Eq. dUi. The intensities of any local magnetic fields (which 
will be assumed imposed in the z-direction) can be then esti- 
mated from {bi}. 

This scheme can be generalized to more complex graphs 
by enlarging the accessible area C as was shown in Ref. 101 ■ 
In that work the Hamiltonian tomography of a general graph 
formed by a network of spin- 1/2 systems was found to be pos- 
sible if C infects the entire graph. The infection process is 
defined as follows. Starting with a subset C of a larger set of 
nodes V, suppose that all nodes in C possess, i.e., are infected 
with some property. This property then spreads and infects 
other nodes according to the following rule: an infected node 
infects a healthy (uninfected) neighbor if and only if the lat- 
ter is the unique healthy neighbor of the former. If eventually 
all nodes are infected by this process, the initial set C is re- 
ferred to as an infecting subset. Figure [T] depicts the infecting 
process with a simple example. We will see below that in 
general, although this requirement of infection is a necessary 
condition of the graph, it is not always a sufficient condition 
for Hamiltonian tomography under arbitrary measurements. 
In particular, we will show that there exist graphs that are not 
amenable to tomography under the spectrally restricted mea- 
surement assumptions employed in the current work. 
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FIG. 1: (Color online) An example of graph infection, (a) Initially, 
three colored nodes in the region C are 'infected', (b) Since the node 
V is the only uninfected node among the neighbors of /i, it becomes 
infected as time evolves, (c) In a similar manner, u' becomes infected 
by . (d) Eventually all nodes will be infected sequentially. 



III. SPECTRALLY RESTRICTED HAMILTONIAN 
TOMOGRAPHY FOR PSEUDO-SPIN NETWORKS 

We now discuss an extension of the gateway scheme for 
Hamihonian tomography of a subsystem, given access to a 
restricted set of spectral measurements and some short time 
subsystem coherence. 

We consider a pseudo-spin network with XY-type interac- 
tions and local external magnetic fields, namely 

= ^ c™„((T+a^ +a„cr+) + ^ fem^"™. (5) 

{m,n)£E meV 

This defines a graph G = (V, E) with V the set of pseudo- 
spin sites and E the links defined by the spin hopping be- 
tween sites. The Hamiltonian parameters are the coupling 
strengths c„i„, and the energy gaps 2fe,„ due to the Zee- 
man terms, ^6,„Z,„. We employ the notation, Xi,Yi, and 
Zi, for the standard Pauli matrices throughout this paper: 
(T^ = (l/2)(X„i + lYm) are thus the raising and lowering 
operators for the m-th pseudo-spin. 

For a network of pigments such as that considered later 
in this paper (Sec. UVl l. the pseudo-spin sites are individual 
molecules with pseudo-spin states |0), |1) corresponding to 
the ground and first excited electronic states, and with energy 
gaps 2b,n, while the links are given by the matrix elements 
of coupling between transition dipole moments on different 
molecules. This corresponds to the usual Heitler-London de- 
scription of excitonic coupling between pigments i29f|. Since 
we restrict the analysis here to short times during which the 
dynamics are coherent, we do not explicitly include other de- 
grees of freedom here (but see discussion in SeclVTi. 

Given a finite window of quantum coherence of subsystem 
dynamics, we may develop a variant of the gateway Hamilto- 
nian tomography scheme via a set of spectral measurements 
at short times. This is possible with the following set of as- 
sumptions. 

1 . The network topology is known. That is, the set of inter- 
acting pairs of sites , which plays a dominant role in the 
overall dynamics, is known without precise information 
on the values of the coupling strengths, c,„„. The latter 
may, without loss of generality, be assumed real. 

2. The sign of each Cmn is known, but not the magnitude. 



3. The energy gaps between the two pseudo-spin levels are 
known for the specific sites that we need to access. 

4. Single site excitation is possible when we have the in- 
formation on energy gaps. 

5. The energy eigenvalues of the system are known. 

6. Measurement in the energy eigenbasis {|£'j)} is possi- 
ble, i.e., the probability of finding an eigenstate \Ej) in 
a state with a single site excitation |n), i.e., |(_Bj|n)p, 
can be measured. 

Assumptions 1 and 2 are the same as in the original scheme 
described in the previous section. Assumption 3 is necessary 
for the single site excitation in Assumption 4. Assumptions 
4 and 5 differ from the assumptions of the original gateway 
scheme of Refs. iH^I], which required waiting for a signal to 
travel back and forth in the chain/network (Sec.|IIli. Since that 
procedure requires long coherence times, in situations where 
the time over which the quantum dynamics are coherent is 
limited, a global measurement in the energy basis provides an 
alternative route to acquire information about the subsystem 
before the excitation decays, provided that such a measure- 
ment may be implemented on a fast enough timescale. With 
use of such a global measurement, the term "to access the site 
ri" then gains a slightly different meaning, namely "to prepare 
a state" or "to excite the molecule" at the site n, rather than 
to "to measure at site n" as was implicitly understood in Ref. 

a. 

The motivation for this measurement in the eigenbasis is 
that all necessary information for the gateway scheme are the 
sets of {Ej} and {Ej\n) for all n ^ C. However, as we 
discuss below, the fact that the quantities obtained from such 
measurements are the modulus of {Ej\n) gives rise to modi- 
fications to the choice of sites that should be accessed and the 
class of graphs to which the scheme is applicable. 

For ID chains, it still suffices to access an end site, as in 
the original proposal in Ref. iQl]. If the graph derived from 
the pseudo-spin network has branches without loops, the end 
sites of all branches should be accessed. Figure |2a) shows 
an example of such a situation. If we set the global phase 
by {Ej\l), then measurements can only give the modulus of 
{Ej\5) (and {Ej\8)), without their relative phases. At site 5 
we have 

{E,-h){E,\5)=c^,{E,\4). (6) 

Summing up the modulus squared of this equation over j, we 
can find c1^. With the assumed knowledge of the sign, we then 
obtain the value of C45, which can then be used to obtain the 
value of |(_Ej|4)|. The procedure is repeated until we reach 
the branching site, i.e., site 3 in Fig. |2ja), where two branches 
meet. The coupling strength between sites 3 and 4 in Fig. 
12 a), for instance, can be obtained by evaluating | • ^ of 
{Ej-bi){E,\4) = C34(-Bj|3)+C45(i;j|4), resulting in = 
Y^jiEj — 64)^|(£'j|4)p, from which 64 can be obtained as 
before. (See text after Eq. (|4]i.) 

If there is a loop in the graph, all sites n that contribute the 
loop need to be accessed in order to determine | (£'j |n) |. The 
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(a) (b) 




FIG. 2: (Color online) Two examples of graphs to which the present 
spectrally restricted Hamiltonian tomography scheme can be applied, 
(a) For graphs with branches, the end pseudo-spin of each branch 
should be accessed, i.e., this pseudo-spin should locally excited be- 
fore the measurement in the |ii'j)-basis. In this example, sites I, 5, 
and 8 (each encircled with red lines) need to be accessed. A mea- 
surement at one of them defines the global phase, e.g., the site 1. The 
coefficients |(i7j|n)| are measured after exciting the red encircled 
sites at the other end of the graph, (b) If there is a loop in the graph, 
then in addition to the end sites of branches emerging from the loop 
(sites 1, 6 and 8), the remaining sites contributing to the loop (e.g., 
site 3, also encircled in red) need to be accessed. Note that in both 
examples we need to access a larger set of sites than the smallest in- 
fecting set. For (a), the smallest infecting set contains only two sites, 
e.g., 1 and 5, and for (b) three sites, such as 1, 3, and 6. 



(a) (b) 




FIG. 3: Graphs containing more than one loop. In both examples 
here, the number of sites is less than the number of edges, rendering 
the Hamiltonian unestimable with the spectrally restricted approach 
using measurements in the energy basis. 




FIG. 4: (Color online) An example of the set of pigment molecules 
to be excited in a small pigment-protein complex. Exciting those 
at sites 1, 5, 6 and 7 (encircled by red dashed lines) individually is 
sufficient to determine all the Hamiltonian parameters. 



necessity of knowing all | (E'j | n) | for the loop-forming sites n 
derives from the requirement of having sufficient equations to 
determine the couplings. If branches extrude from the loop, 
the access sites should be chosen to be the end sites of these 
branches, just as in the case of simple graphs having branches. 
(See Fig. Eb).) 

These examples show that the present variant of the gate- 
way scheme for Hamiltonian tomography cannot be applied 
once there are two or more loops in a connected graph. The 
constraint on the available measurement given by Assumption 
5 above poses a further condition on the network topologies to 
which the current scheme is applicable. We illustrate this for 
two simple examples of graphs in Fig. [3] Because the original 
gateway scheme is based on the fact that only one unknown 
term, e.g., Cmn{Ej\n), appears in the equation deriving from 
the factor {Ej\H\m), we are able to obtain the value of one 
new coupling strength, c,„„, by using the previously obtained 
knowledge on parameters for sites other than n. The property 
of infection then guarantees that all coupling strengths can be 
estimated recursively in this manner, provided that all coeffi- 
cients {Ej\m) are known for all sites m G (7, as well as all 
eigenvalues {Ej}. However, when there is no information on 
the relative phase of {Ej |m), the number of sites needs to be 
larger than or equal to the number of edges. The graphs in Fig. 
[3]do not fulfill this condition and as a result their Hamiltonians 
cannot be estimated with the current approach. 



IV. APPLICATION TO A MOLECULAR NETWORK 

We now apply the gateway scheme with restricted spec- 
tral access outlined above for a network of pseudo-spins to 
Hamiltonian tomography of a network graph with nontrivial 
topology mimicking a pigment-protein complex. As an exam- 



ple, we take the geometry of the seven coupled pigments em- 
bedded in the Fenna-Mathews-Olson protein luTT and use the 
dominant electronic couplings between pigments as in Refs. 
ifsojlslll . This leads to the network graph shown in Fig. H) 

Recent investigations have demonstrated that electronic 
quantum coherence in such a pigment-protein complex per- 
sists for several hundreds of femtoseconds even at physio- 
logical temperatures ll23l - l26ll . This means that application of 
the present Hamiltonian tomography scheme is restricted to 
measurements on the timescale of a few hundred femtosec- 
onds. We note that the presence of quantum coherence does 
not necessarily imply purely unitary dynamics. If the evolu- 
tion is indeed unitary, the values of |(i?j|n)| are constant in 
time. However, in the presence of dissipation, the measured 
values of |(i?j|n)| may vary. Provided that the measurement 
can be made within the timescale in which the dynamics of 
\Ej) may be characterized by a phenomenological factor Tj, 
then this time dependence would be reflected in measurement 
of time dependent coefficients \{Ej\n)\cxp{—rjt/2). Mea- 
suring these quantities at various times within the relevant 
timescale would then allow estimations of the values | {Ej\n) \ 
by extrapolation back to i = 0. CuiTent technology allows 
controlled shaping of pulses with time duration 10-20 fem- 
toseconds, suggesting that such an estimation might be feasi- 
ble. 

Using the data of {Ej} and {{Ej\ni) (m e C)} and fol- 
lowing the procedure described above with Eqs. ©-(U), we 
can then construct the matrix elements of the symmetric ma- 
trix corresponding to the Hamiltonian in the one-excitation 
subspace. 

Since the network in Fig. |4]contains a loop formed by four 
sites, we need to access at least four sites, 1, 5, 6, and 7, as 
described in Sec |III1 Let us follow the estimation procedure 
again briefly for clarity. Suppose that we start from that 
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is, we set the global phase of \Ej) so that {Ej\l) are real and 
positive for all j. Then, using Eq. (|2]i, we have 

Y}{E, -h,){E,\l)\^ = c,2{E,\2), (7) 
j 

where hi can be known from Ej and {Ej 1 1) , thus the left-hand 
side of Eq. Q is equal to {Ej 1 1)^ — b\. The estimation 

process then proceeds to site 4 according to Eq. (|3]l, obtaining 
the values of C34 and {Ej\4) with a correct phase. With the 
measured values of KSjiSj)!, KGlSj)], and |(7|£'j)|, we can 
then make use of the following set of equations to obtain the 
coupling strengths: 

{E, - hi){E,\A) - C34(i^,|3) = C45(i?,|5) + c„{E,\7), 
{E, - h){E,\b) = c,,{Ej\4) + c,e{E,\6), 
{E, - h){E,\Q) = C56(i?,|5) + C67(i5,|7), 
{E, - bj){E,\7) ^ C67(i?,|6) + c„{E,\4). 

(8) 

Summing up the modulus squared of each equation over j 
gives four equations with four unknown parameters, c^g, Cgg, 
c^y, and €47. Together with the a priori knowledge on the 
signs of the set {c„} that was assumed initially, all coupling 
strengths Cij can now be estimated. The remaining parame- 
ters, i.e., 65, h^, and 67, may be evaluated from | {Ej |n) | (n ~ 
5,6,7). For example, &5 = (5|H|5) = Ej\{Ej\5)\'' . 
Thus all parameters of the Hamiltonian have now been iden- 
tified, despite the lack of the precise information about the 
phase of {Ej |n). 

V. DISCUSSION 



of such a spectrally restricted measurement, we then showed 
that by choosing the right set of accessible (i.e., spectroscop- 
ically excitable) sites, the Hamiltonian of a given network of 
pseudo-spins can be estimated. These constraints on measur- 
able quantities, in particular the lack of feasibility of mea- 
surements in the site basis, were found to modify the require- 
ments for the graph properties. While the Hamiltonians of 
one-dimensional chains are still estimable by accessing, i.e., 
preparing a state at the end site only, as in the original gateway 
scheme, we now find that, in general, the set of accessible sites 
needs to be larger than an infecting set. Furthermore, there ex- 
ist networks to which the current estimation scheme cannot be 
applied, regardless of the choice of accessible sites, because of 
the detailed structure of the network topology. This major dif- 
ference results from the constraint on feasible measurements 
imposed here and raises interesting questions for the interplay 
between network topology and measurement capabilities in 
Hamiltonian identification schemes in general. 

Application of this spectrally restricted Hamiltonian to- 
mography approach to a small scale network of molecular 
pigments indicated that provided the spectral measurements 
can be made on a timescale significantly shorter than the 
characteristic time for loss of coherence, the Hamiltonian to- 
mography approach may yield useful estimates for parame- 
ters of the electronic Hamiltonian describing excitonic energy 
transfer through a well-characterized aggregate of pigments or 
pigment-protein complexes. We note nevertheless that realis- 
tic application to such systems will require extension of the 
current approach to include dissipation and decoherence in a 
more quantitative manner, e.g., as in Refs. 14, J., 12, 15, 32j- 
[34I1 . to account for the lack of unitarity that is associated with 
the subsystem dynamics despite the appearance of quantum 
coherences at short times. 



We have developed an extension of the gateway scheme of 
Hamiltonian tomography to estimation of Hamiltonian param- 
eters for subsystems of complex quantum systems that show 
coherent dynamics for a limited period of time. We circum- 
vent the problem of decay of quantum coherence preventing 
the observation of reflections of the injected signal that was re- 
quired in the original scheme of Ref. jlj] by employing instead 
a measurement in the {|i?j)}-basis. Assuming the feasibility 
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